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I. INTRODUCTION 



Since the pioneering work in Ref. the rich phase structure in the lattice Wilson fermion 
has been extensively studied As is well-known [3], Wilson fermions bypass the no-go 

theorem and produce a single fermionic degree of freedom by breaking the chiral symmetry 
explicitly. This leads to an additive mass renormalization and requires fine-tuning of a mass 
parameter for a chiral limit. Furthermore at finite lattice spacing, there emerges a parity- 
broken phase (Aoki phase) The full phase diagram reflects the masses possessed by each 
of the original doublers. As seen from this fact, the main reason for the emergence of the 
parity-broken phase is that the Wilson term gives a species(taste)-sensitive mass to produce 
a mass splitting of species as well as breaking the chiral symmetry. The understanding of 
the parity-broken phase structure is not only useful for simulations with Wilson fermions, 
nit also gives practical information for the application of overlap jo], [lo| and domain-wall 
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12| fermions, both of which are built on the Wilson fermion kernel. Indeed it is shown in 



13| that the domain- wall fermion also possesses a complicated parity broken phase diagram 
for a finite size of the extra dimension. 

On the other hand, no parity-broken phase structure is observed in staggered fermions 
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16| with their exact chiral symmetry. However things could be changed if we introduce 



a taste-sensitive mass term, which we refer to as a taste-splitting or flavored mass in this 
paper. Adams recently established theoretical foundation of the index theorem with stag- 
gered fermions 17l| a nd presented a new version of the overlap fermion constructed from the 
staggered kernel 18|, ll9| . He introduced a taste-splitting mass term for the spectral flow to 
detect the index correctly. This mass term assigns positive and negative masses to tastes 



depending on their flavor-chiralities. After these works the present authors 
defined the index in the naive and minimally doubled fermions 
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20| successfully 



24j and presented new 



versions of overlap fermions by implementing the flavored mass terms [25|. It is natural to 
consider the phase diagram for these fermions with the mass splitting of the tastes since it 
is also useful for the practical application of their overlap versions as well as themselves. 
In this paper we study the parity-broken phase structure for naive and staggered fermions 



with t 



le flavored mass terms. We use the two-dimensional lattice Gross-Neveu models 



27| as toy models of QCD. We develop the generalized staggered Gross-Neveu model 



including two types of four-point interactions to study the staggered phase structure. We 
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solve the gap equations for the large limit and obtain the phase boundaries in the M-g"^ 
plane. We show the Aoki phase exists both in staggered and naive cases reflecting the mass 
splitting of tastes. In the naive cases there are varieties of the phase diagram depending 
on linear combinations of two types of the flavored masses. This elucidation of the phase 
structure can contribute to the practical application of these fermions and their overlap 
versions. We also discuss the continuum limits of these Gross-Neveu models. We show that 
we can take the chiral continuum limit with the associated number of massless fermions 
by flne-tuning a mass parameter and two coupling constants. It indicates that, in Lattice 
QCD with the staggered fermions with the Adams-type 18| or Hoelbling-type 19|] flavored 
masses, we can obtain the two- or one-flavor massless fermions in the chiral and continuum 
limit by tuning the mass parameter. They can be less expensive than Wilson fermions in 
lattice QCD simulations. 

In Sec. [TTl we study the parity broken phase diagram by using the naive Gross-Neveu 
model with the flavored mass. In Sec. Illll we propose the generalized staggered Gross-Neveu 
model and study the phase diagram. In Sec. IIVI we investigate the continuum limit of these 
models and discuss the flrst order phase phase boundaries in the phase diagram. Section |V] 
is devoted to a summary and discussion. 



II. NAIVE GROSS-NEVEU MODEL 



In this section we investigate the phase diagram for naive lattice fermions with flavored 
mass terms by using the d = 2 Gross-Neveu model, which has lots of common features with 
QCD. Let us begin with the lattice Gross-Neveu model with the flavored mass term, which 
is given by 

1 - g"^ - - 

n,fi n 

+ J2MMSnm + {Mf)^Jij^, (1) 
n,m 

where /i stands for fi = 1,2, n = [rii, 77,2) are the two dimensional coordinates and ipn stands 
for a A^-component Dirac fermion fleld {ipn)j{j = 1,2,..., A^). We note the bilinear ipip means 
Ylf=i'^j'^j- 9^ corresponds to the 't Hooft coupling. M is a usual mass assigning the same 
mass to species while {Mf)n,m is a flavored mass assigning different masses to them. Here 
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we define the two dimensional gamma matrices as 7i = 72 = (J2 and 75 = (T3. We make 
all the quantities dimensionless in this equation. By introducing auxiliary bosonic fields (T.„, 
TTn we remove the four-point interactions as 

n,ii n,m 

+ ^ ^{{<^n - Mf +7^1) + Man + il5Trn)M (2) 

n n 

By solving the equations of motion, we show the following relation between these auxiliary 
fields and the bilinears of the fermion fields 

an = M - ^^ij, (3) 
TTn = -^'^ilbi^- (4) 

These relations indicate how a and tt stand for the scalar and pseudo-scalar mesons. After 
integrating the fermion fields, the partition function and the effective action with these 
auxiliary fields are given by 

Z = I n(rf^nCivr„)e-^^-('^'^), (5) 

n 

Scs{an, 7r„) = ^{{an - Mf + vr^) - Tr logZ^^^^, (6) 
9 



with 
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Dn,m = {an + il5'^n)Sn.m + "^(^n+Z^,™ - ^n-^.m) + {Mf)n,m- (7) 



Here Tr stands for the trace both for the position and spinor spaces. As is well-known, the 
partition function in the Gross-Neveu model is given by the saddle point of this effective 
action in the large limit. We denote as 5"„, 7r„ solutions satisfying the saddle-point 
conditions 

an , T^n an, T^n) „ /„\ 

ban 

Then the partition function is given by 

Z = e-^^«(^'*). (9) 

By assuming the translational invariance we define the position-independent solutions as 
ao = o"o and ttq = ttq Then we can factorize a volume factor V = J2n ^ effective action 
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as 



S'cff = \^S'cflf(cro,7ro), (10) 
4ff(ao,vro) = ^{{cro-Mf + 4)-^TT\ogD. (11) 

We can write Tr log D in a simple form by the Fourier transformation to momentum space 

/d'^k \ 
log[det((To + nsVTo + Mf{k) + i22 7m sin k^)] 

= (0^°S[^^° + + + ^^^^ 

with det being the determinant in the spinor space and = ^^sin^/c^. Mf{k) is the 
flavored mass represented in momentum space. Now the saddle-point equations are written 
as 

SS,s_ {(To-M) f d^k ao + Mf{k) 

5^0 J (27r)2(ao + M;(A;))2 + vr2 + s2 ' ^ ^ 

SS^ ^ no _ f jfk vTo ^ 

^vro J i27ryiao + Mfik)y + 71^ + 3^ ' ^ ^ 

In this section we consider two types of the flavored mass for the naive fermion 

Mf{k) = cosA:icosfc2, (15) 
Mf\k) = -(cos /ci + cos A;2)(l + cos /ci cos /i;2)- (16) 

Such mass terms were flrst introduced in the minimally doubled fermion by using the point- 



splitting method 251. Then these were introduced also for the naive fermion to consider 
the index theorem and a new type of overlap fermions [20(. Studying the phase diagram 
with these flavored mass terms not only contributes to understanding the overlap versions 
but also helps to understand the staggered case in the next section. Here ctq and ttq are 
determined as cro(M, (7^), 7ro(M, gf^) from the saddle-point equations once the values of M 
and are flxed. 

Let us look into the phase structure with respect to parity symmetry. The order parameter 
of this symmetry is ttq, which can take zero or non-zero values depending on values of M 
and g^. Parity symmetry is spontaneously broken for the non-zero cases ttq 7^ 0. The phase 
boundary is determined by imposing ttq = on Eq. ( fT3l)(fT^ after the overall ttq being 
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removed in Eq. (fT^ . Then the conditions for the phase boundary, so-called gap equations, 
are given by 

J (27r)2(ao + M^(A;))2 + 52' ^ ') 

^ f d'k 1 

J (27r)2(ao + M;(A;))2 + s2' ^ ^ 

with Mc being the critical value of M. As we will check later, this phase boundary is a 
second-order critical line. Here we derive the parity phase boundary Mc{g'^) as a function 
of the coupling g'^ by getting rid of the chiral condensate ctq from these equations. We will 
calculate the parity phase boundaries for three cases of the flavored masses Mf \ Mf and 
Mf^ +Mf. 



A. Mf 

The lattice fermion action with this flavored mass assigns the positive mass m = 1 to 
two species with the momentum (0, 0)(7r, vr) and the negative mass m = — 1 to the other two 
species with (0, 7r)(7r, 0). Before calculating Mc{g'^) numerically, we can anticipate the phase 
structure from the symmetry of the gap equations. To see this we replace ki hj tt — ki in 
(US]) and (IH]) for M^^\ Then the equat ions are converted into 



ao + M f <Pk -ap + Mf\k) 



TTq f d?k TTo 



(19) 
(20) 



g^ J {2Tif {^^ao + Mf\k)f + Til + s^' 

Thus, if ((To, VTo) are solutions for (M, g"^), (— (TQ) ttq) are solutions for (— M, g"^). It also 
means, if (Mc, g"^) is a critical point, (— M^, g"^) too. We can anticipate the phase diagram 
for this case is symmetric about M = 0. Now we derive the parity phase boundary Mc{g'^) 
numerically for M^p (k) = cos ki cos k2- The phase diagram for this case is depicted in Fig. [H 
A stands for the parity symmetric phase ttq = and B for Aoki phase vto 7^ 0. In the large 
coupling region there are two phase boundaries while there are four phase boundaries in the 
weak coupling region. The left and right cusps correspond to two species (0,0)(7r,7r) with 
the positive mass (m = 1) and the other two (0,7r)(7r,0) with the negative mass (m = — 1) 
respectively. It reflects the mass splitting of species given by the flavored mass Mf \ Here 
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FIG. 1: Aoki phase structure for the naive fermion with the flavored mass My' . The left and right 
cusps are related to two species (0, 0)(7r, vr) with m = 1 and the other two (0, 7r)(7r, 0) with m = — 1 
respectively. A and B stands for parity-symmetric and -broken phases. 

we note we obtain the same result for — ^j-^^ except that the species (0,0)(7r,7r) live at the 
right cusp and the other two live at the left. It means the sign of the this flavored mass is 
irrelevant for the spectrum of the Dirac operator or the associated Aoki phase. 



The lattice fermion action with this flavored mass assigns the positive mass (m = 2) to 
one of four species with the momentum (0,0), zero mass to (0,7r)(7r,0) and the negative 
mass (m = —2) to (7r,7r). To look at the symmetry of the gap equations we replace A;^ by 



Thus, if ((To, TTo) are solutions for (M, g^), {—cto, ttq) are solutions for (— M, g'^). It also 
means, if (Mc, g"^) is a critical point, (— Mc, g"^) too. We can anticipate the phase diagram 
for this case is again symmetric about M = 0. Now we derive the parity phase boundary 
Mc{g'^) numerically for Mj-^''(A;) = (cos /ci -|- cos A;2)(l + cos /ci cos A;2)/2. In the large coupling 
region there are two phase boundaries while there are six phase boundaries in the weak 
coupling region. The three cusps correspond to one of four species (0, 0) with m = 2, two of 



B. m: 




(21) 



(22) 
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FIG. 2: Aoki phase structure for the naive fermion with the flavored mass Mj . The three cusps 
correspond to (0, 0) with m = 2, (0, vr)(7r, 0) with m = and (tt, tt) with m = —2 respectively from 
the left. 

them (0,7r)(7r,0) with m = and the other one (tt, tt) with m = —2 respectively from the 

(2) 

left. It reflects the mass splitting of species given by the flavored mass '. 

c. uf^ + uf^ 

The fermion action with this flavored mass assigns the positive mass m = 3 to one of 
species with the momentum (0, 0) and the negative mass m = — 1 to the other three species 
with (0, 7r)(7r, 0)(7r, tt). Here we cannot find any relevant symmetry in the gap equations. 
Thus we can anticipate the phase diagram for this case is not symmetric. Now we calculate 
Mc((7^) numerically for M^p + {k) = cos ki cos k2 + (cos ki + cos k2){l + cos ki cos k2)/2. 
The result of the phase diagram is depicted in Fig. [3l It is obvious that it is not symmetric 
about M = 0. In the large coupling region there are two phase boundaries while there are 
four phase boundaries in the weak coupling region. The left and right cusps correspond 
to one of species (0,0) with m = 3 and the other three (0, 7r)(7r, 0)(7r, tt) with m = — 1 
respectively. It reflects the mass splitting of species given by the flavored mass 
Now we can easily modify the phase diagram by choosing the linear combination of Mj^^ 
and Mf. 

We expect these results are qualitatively similar to the phase diagram of the d = 4 
fermion actions with the Non-abelian gauge field like QCD except for the number of species 
associated with each cusp. In the end of this section we check the mass of the vr-meson 
becomes zero on the critical line Mc{g^). As is well-known, the correlation length gets 
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FIG. 3: Aoki phase structure for the naive fermion with the flavored mass My' + My> . The left 
and right cusps correspond to (0,0) with m = 3 and (0, 7r)(7r, 0)(7r, vr) with m = —1 respectively. 

infinitely large in the vicinity of the second and higher phase boundaries, which leads to 
massless dynamical degrees of freedom. In the case of lattice QCD with chiral-symmetry- 
broken fermions like Wilson fermion, the fine-tuning of the mass parameter to the 2nd order 
phase boundary leads to the chiral limit with massless quarks and massless pions regarded 
as Goldstone bosons due to the spontaneous chiral symmetry breaking. Thus it is quite 
important to verify it. We can show the mass of vr^ becomes zero on the phase boundaries 
as 



The zero mass of the pion means the phase boundary we derived is the second-order critical 
line. We can also check the order of the phase boundaries by depicting the potential for (Tq 
and ttq as we will discuss in Sec. IIVI 



In this section we investigate the phase diagram for staggered fermions with the Adams- 
type flavored mass term by using the d = 2 Gross- Neveu model. To study the parity broken 




0. 



(23) 



III. STAGGERED GROSS-NEVEU MODEL 
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phase structure we propose the generahzed staggered Gross-Neveu model with the 75-type 
4-point interaction, which is given by 



n,fj, n 

- [(E^2^+^^2Ar+^)' + (E^(-l)^''"^'^2A^+^^2Ar+A)'], (24) 

Af A A 

where we define two-dimensional coordinates as = 2N'+A with the sublattice A = {Ai, A2) 
(^1,2 = 0, 1). Xn is a one-component fermionic field. (—1)^1+^2 corresponds to the natural 
definition of 75 for this fermion which is expressed as T^^ = 75 ® 75 in the spinor-taste 
expression. 77^ = (— corresponds to 7^. As the fiavored mass term we choose 
the Adams-type one, which is given by 

My = r5r55~l®75 + 0(a) (25) 

with the following chirality matrix Fs 

= -tviV2j2CiC2, (26) 

sym 

C, = l{T, + T^,) (27) 

where is the usual translation operator. (The chirality matrix in general dimensions is 
defined as = —{iY^'^rii ■ ■ ■ ?7rfXlsym^i ' ' " ^d-) This mass term assigns the positive mass 
(m = +1) to one taste and the negative mass (m = —1) to the other depending on ± 
eigenvalues for FsFss which we call the fiavor- chirality. With bosonic auxiliary fields aj^f, 
TTj^j-, the action is rewritten as 



5" = ^ ^Vf^XniXn+f, - Xn-^) + ^ XnMfXr. 



2 



+ ^ - M)2 + vr^) + ^ X2X+a{ctm + ^(-l)^^+^^vr^)X2Ar+A, (28) 

After integrating the fermion field, the partition function and the effective action with these 
auxiliary fields are given by 

Z = j (I?aAr2^7r^^)e"^^^fl('^'"\ (29) 
-^eff = ^ + - log D, (30) 

AT 
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with 



\Ai+A2 



(5„, 



6^ 



(31) 



The process from ([8]) to ( ITT]) in the case of the naive fermion is common with this staggered 
case. We again denote as (Tq and ttq the position-independent solutions of the saddle-point 
equations. In this case, however, it is not straightforward to derive the Trlog-D with the 
Dirac operator fl3T]) in the effective action Eq. f|TT]) . In order to estimate this trace logarithm 
we first obtain the determinant of the Dirac operator in the sublattice space, which means 
the determinant in the spinor and taste spaces. Here we express the sublattice structure as 
a multiplet field XN composed of the four one-component fields as 



Xiii 

Vxiv/ 

where we mean i = 2M, ii = 2J\f + (1, 0), iii = 2A/' + (0, 1) and iv = 2A/' ■ 
estimate the trace term 



(32) 



Tr logD = V 



hgdet{{D{k))ab, 



:i, 1). Now let us 



(33) 



where a, b stand for the index of the four sublattices running from i to iv. Here det means 
the determinant with respect to the sublattice. The Dirac operator is given by 



{D{k))ab = 0-oSab + 



+ 



ZTTo 



+ i 



+ 



ki k2 
cos — cos — 
2 2 



/ 



V 



I sm 



i sin ^ 
i sin f 



I sm 



k2 



-z sin f 
z sin ^ 



-i sin Y i sin ^ 



/ 



(34) 
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Then detD is given by 

detiD{k))ab = {al + Til + s^f - 2clclial - 4 - s^) + 44 

= ((ao + cica)^ + TT^ + s^){{ao - CiCa)^ + tt^ + s^), (35) 

where = sinA;^/2, = Yln^fi^ '^m ~ cos/c^/2. It is notable that this determinant is 
expressed by the product of the two determinants of the naive fermions with the flavored 
mass ±Mf{kf,/2). Now we can explicitly write the saddle-point conditions satisfied by ao 
and ttq as 

(To - M _ ^ f dk'^ ao{a^ + tt^ + s^) - c\c\aQ ^^^^ 



J (27r)2 ((ao + CiCs)^ + vr^ + s^){{ao - CiCs)^ + tt^ + s 
1 r dk'^ 0-^ + 5^ + c?c2 



J (27r)2((ao + CiC2)2 + 7r2 + s2)((ao-CiC2)2 + 7r2 + s2)' 

^ f ^o(^0 + + S^) + cfclTTo 

J (27r)2((ao + CiC2)2 + 7r2 + s2)((o-o-CiC2)2 + 7r2 + s2)- ^ ^ 

By multiplying —1 to the first equation, we see (— O"0) ^0) are solutions for (— M, g"^) if (ctq, 
ttq) are solutions for (M, (7^). It also means, if (M^, (7^) is a critical point, (— M^, (7^) too. 
The phase diagram will be symmetric about M = 0. The parity phase boundary Mc{g'^) in 
this case is derived by imposing ttq = in (l36l) ( 1371) after the overall ttq being removed in the 
second one. Then the gap equations are given by 

Mc , f de 2cl4ao 

'1^2J -r 'iQ -r o 

_ -r -r ^1^2 , , 

V (27r)2((cro + CiC2)2 + 7r2 + s2)((ao-CiC2)2 + 7r2 + 52)- ^ ^ 

By removing ao in these equations, we derive the phase boundary Mc^g"^). The result is 
shown in Fig. HI 

Here again A stands for the parity symmetric phase (vto = 0) and B for Aoki phase 
(tto 7^ 0). In the large coupling region there are two phase boundaries while there are four 
phase boundaries in the weak coupling region. The left cusp corresponds to one of two tastes 
with m = 1, and the right corresponds to the other taste with m = — 1. Thus the phase 
diagram reflects the mass splitting of tastes given by the Adams-type flavored mass. We 
also check the pion mass becomes zero on the second order phase boundary as 

°^ — )|m=m, = V—^\m=m, = 0. (40) 

Now let us consider the parity phase structure in the d = A QCD with the staggered fermion 
with this flavored mass. Considering the case of the Wilson fermion we can speculate it 
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FIG. 4: Aoki phase structure for the staggered fermion with the Adams-type flavored mass FsFss. 
The left and right cusps correspond to one of two tastes with m = 1 and the other with m = — 1. 
A stands for a parity symmetric phase and B for Aoki phase. 

is qualitatively similar to our result for the d = 2 Gross-Neveu model except the number 
of species associated with each cusp. In the four dimension, four tastes in the staggered 
fermion with the Adams-type flavored mass split into two with positive mass and the other 
two with negative mass depending on their flavor- chiralities. Thus each of the cusps in the 
phase diagram will correspond to two tastes. If we consider another type of the flavored 



mass term proposed in 



19|, the four tastes are split into one with positive mass, two with 



zero mass and the other with negative mass. If we can take the chiral and continuum limit 
around the cusps, we obtain the two- or one-flavor staggered fermions with tuning only the 
mass parameter, which will be numerically faster than Wilson fermion. Thus the question 
here is whether we can take the massless continuum limit. We will discuss this point in the 
next section with starting with the case of the naive fermion. 

IV. CONTINUUM LIMIT 



In this section we discuss the continuum limit of the naive and staggered Gross-Neveu 
models with the flavored masses discussed in Sec. [TT] and Sec. IIIII This analysis gives us 
important informations on the continuum limit of the d = A QCD with these fermions. As 
is well-known, the chiral symmetry is realized in the effective potential of the Gross-Neveu 
model as the 0(2) rotational symmetry about ao and ttq. The purpose here is to figure out 
the fine-tuned values of the mass and couplings to recover this symmetry for a — 0. We 
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note in order to take the chiral and continuum limit in this model, we need to introduce 
u™ independe. coupUn.. A a„d A a. we wiU ,ate. The *a.e.y . .o expand .Ke 
fermion determinant in the effective potential with respect to the lattice spacing a following 



the process in 



3|. 



We first consider the case of the naive fermion with one of the flavored masses M^^ = 
cos fci cos /c2- The effective potential in this case with the lattice spacing being explicit is 
given by 

^ = / -j-^log[(go + -cosfciQcosfcaa)^ + + X] 2^^" ]- (^^) 
Now we divide the terms in the determinant / into (9(l/a^) and 0{l/a) parts as 



I {Do, D,) = ^ \og[Do + Di], (43) 

^ sin^ /c„a , a,, , /a + cos /cia cos /i;2a\ 2 

Do = 22 ^ + (^0 - + 4 + • 44 

a, /a + COS fcio COS A;20\ /.^x 

= 2 fTo - - . 45 

a \ a J 

where we introduce a constant a since there is arbitrariness about how to divide the terms 
into O{l/o?) and 0{l/a) parts. This is determined by which cusp you choose in Fig. [1], or 
equivalently which species you want to make massless in the continuum limit. Here we fix 
a = —1 which is related to the left cusp or the continuum limit with the massless species 
(0,0) and (7r,7r). (With a = 1 we can discuss the other cusp while we will discuss a = 
in the end of this section.) Here we use the shifted definition of ctq as ctq + 1 /a- — > ctq for 
simplicity for a while. Then the effective potential with this shift is given by 

^ v-^ sin^ fcuO 9 9 / — 1 + cosA;iacosA;9a\ 2 

/-I + cosfciacosfcsax 
Di = 2ao j . (48) 
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We expand / by Di/Dq or equivalently by the lattice spacing a, 

/ = /o + ^/„, (49) 



n=l 



^0= I 77^ log 1^0, (50) 



n y_/,(27r)2Do" 


(ri>l), 


V / \nn—2 

n 






(-1 + COS^i COS ^2)" 


(2^)ME,sin2e. + ( 


;-l + cos^i cos6)2 + a^{a^ + tt^))'' 



(51) 

where we introduce the dimensionless momentum C,fi = k^a. For a — )■ 0, only the Jq, /i and 
I2 remains nonzero. /o(a — 0), /i(a — )■ 0) and 12(0 — 0) are given by 

/o(a -> 0) = C,{al + vTo^) - -(a^ + 4) log + [C^ = 0.367), (52) 

ZTT e 

/^(a ^ 0) = (Ci = -0.446), (53) 

hia ^ 0) = -2alC2 (C2 = 0.201). (54) 

From here we basically do not care about the 0{a) corrections. Here we show the explicit 
values of Cq, Ci and C2 since they will be essential for the discussion later. The details of 
the calculations are shown in Appendix I A 1[ Now let us discuss the continuum limit of this 
theory. Including all the nonzero contributions for a — )■ 0, the effective potential is given by 



M + l/a , 2^^^ , ^ 1 ^ , 1 ,_^2\„2 

2 I _2 



ScS = -{ -2 + -CAao+ (—^ - Co + —\oga')7r. 



+ - C^o + 2C2 + log a') + ^{al + 4) log (55) 

This indicates we need two independent couplings g^, to recover the 0(2) symmetry 
toward the continuum limit. In addition, getting rid of the ctq linear term leads to the 
massless limit. Then the natural fine-tuned parameters for the chirally symmetric continuum 
limit without 0{a) corrections are given by 

M = — ^Ci - 1, (56) 



4C2gf+l' 



9l = 777^, (57) 
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where Eq. (156|) is obtained by imposing the coefficient of (Tq and Eq. (1571) is given by imposing 
the coefficients of Uq and tTq coincide. To consider a renormahzed theory with the chiral 
symmetry we introduce the scale parameter (A-parameter) as 



Aa = exp 



TT 

ttCo - 27rC2 - 7-^ 



(58) 



With the natural fine-tuning (1571) . this definition of A leads to the coupling renormalization 
including a given by 

_L = c„-2C. + J^log(J^). (59) 

Here we need to keep A finite when we take the continuum limit a — > 0. Then the renor- 
mahzed effective potential with the chiral symmetry in the continuum limit is given by 



^e. = ^K^ + vro^)log^^ (61) 



We note the fine-tuned point {M{gl), gl{gl)) in ( l56|)(l57|) specifies the line along which the 
continuum limit should be taken. At the minimum of this potential o"o has a nonzero value, 
which corresponds to the spontaneous chiral symmetry breaking. 

Let us look at these fine-tuned parameters in terms of the phase diagram. By this we 
can verify our fine-tuning yields the chiral-symmetric continuum theory. We ffist consider 
the non-zero value of g"^ as g^ = 0.6 to reveal properties of the phase diagram. By hiding 
the lattice parameter with a = 1 the fine-tuned point (M(0. 6), (7^(0.6)) is given by 

M{gl = 0.6) = -0.464, (62) 
gligl = 0.6) = 0.404. (63) 



Now we consider the M-g^ phase diagram with g^ = 0.6. According to the case of the Wilson 
Gross- Neveu model js], the phase boundary has a self-crossing point and the fine-tuned point 
is located slightly inside and below the self-crossing point in the parity symmetric phase. 
Besides the phase boundary naively derived from the gap equations no longer describes 
the true one near the self-crossing point, and we need study the effective potential to find 
the true critical lines including the 1st order ones. Here we will show these situations are 
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common with our cases. The gap equations for the two couphngs are given by 

= aJl - ^] - 2gl ——- ^j-^ , 64) 

Here we come back to the unshifted definition of ctq. In Fig. |5] and Fig. [6] we depict the 
Mc{gl) phase boundary derived from the gap equations ( I64|) (l65l) for gl = 0.6. The latter is 
an expanded one near the self-crossing point with the true phase boundaries. In the both 
figures a crosspoint stands for the fine-tuned point without 0{a) corrections (162|) ( 163|) . It is 
located slightly to the right and below the self-crossing point near the second order phase 
boundary. We note this region is the parity-unbroken phase. The qualitative properties of 
this phase diagram remain toward g^ ^ where the whole structure moves down to g^ = 
with the Ist-order boundaries disappearing. Here the fine-tuned point fl56|)fl57j) gets close 
to the endpoint of the 2nd-order phase boundary at (M, (7^) — )■ (—1,0), which corresponds 
to two species (0, 0)(7r, tt). Thus the continuum limit along this fine-tuned point yields the 
theory with chiral symmetry and two massless fermions, which leads to massless pions as 
Goldstone bosons. 

Now we discuss the first order phase transition. Although it is not essential for our 
purpose because in the limit g'^ the Ist-order phase boundary disappears and the 
entire phase boundary becomes of 2nd order, we can reveal other aspects of our fermions by 
investigating it. As shown in 6| there are two kinds of the 1st order phase boundaries in the 
case of Wilson fermion. One is the parity phase boundary, across which ttq at the minimum 
of the effective potential changes from zero to nonzero. The other is related to ctq, across 
which the sign of ctq at the minimum of the potential changes discontinuously. Now we will 
show both of them exist also in our case. We numerically calculate the effective potential in 
Eq. (142!) and search the minimum of the potential. In Fig. [6] we depict the appearance of the 
1st order phase boundaries. Here we note the true parity phase boundary of 2nd order as 
a blue solid line coincides with the naively derived phase boundary as a blue dotted line at 
the both sides of the self-crossing. Then the 2nd-order one coming from the left converts to 
the Ist-order at some point, which is spilled out from the naively derived boundary. It ends 
at the point encountering the naively derived one again. The Ist-order phase boundary for 
ctq starts from this point, going down straight, and ends at gl = 0. In Fig. [7] we depict the 
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FIG. 5: The naively derived phase boundary M{g'^) for the naive fermion with with = 0.6. 
The fine-tuned point (—0.464, 0.404) as a crosspoint is located near the self-crossing point. 

order parameter ttq as a function of M for some fixed values of around which the order 
changes in Fig. O Here we verify the order of the transition changes from the 2nd to the 1st 
about the point. In Fig. [H]we depict the ctq potential for several values of M crossing the 
(To phase boundary. (Here we can take ttq = since it is the parity symmetric phase.) The 
value of (To at the minimum changes from ctq > —1 io < —1 m. & form of the Ist-order 
phase transition. Indeed the potential describing these Ist-order transitions is also obtained 
by taking account of 0{a) corrections. The contribution from the correction 5Scs is given 
by 

55efr = -^Cal + 2ao{al + 4) [c^ + ^ log ^i±ll^ , (66) 

with C3 = —0.0923 and Ci = —0.0741. We can qualitatively reproduce the above results 
from the effective potential with these corrections. We can obtain the same but reversed 
phase structure for the right cusp by choosing a = 1/a in fH^fHSj) . We also note the 
sign of (To continuously changes at M = 0. It is related with the discrete chiral symmetry 
((To —7- —(To) of the effective action ( l42i) for M = up to a irrelevant sign. This symmetry 
indicates interesting possibility of another continuum limit corresponding to the case of 
Q? = in (l44l) (145!) . We will discuss details on this topic in the end of this section. 

In Fig. [9] and [10] we depict the corresponding figures for the flavored mass 
We take g"^ = 1.2 to make the structure enhanced, where the fine-tuned point for the left 
cusp is given by {M,gl) = (—2.205,0.720). The results are qualitatively the same as the 
previous case. In this case the continuum limit along with the fine-tuned point leads to the 
single-flavor theory with one of the species at (0, 0). 
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FIG. 6: An expanded version of Fig.[5l A blue dotted curve is the naively derived phase boundary. 
The true phase boundaries are composed of the three parts. The fine-tuned point as a cross point 
is located slightly to the right and below the self-crossing point. 
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FIG. 7: The order parameter vro as a function of M for = 0.41, 0.42, 0.43 where the order of 
transition changes from 1st to 2nd in Fig. [H 
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FIG. 8: The cq potential for several values of M crossing the do boundary in Fig. [6l The value of 
(To at the minimum changes from (Tq > — 1 to cJo < — 1 in a form of the Ist-order transition. 
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FIG. 9: The naively derived phase boundary Mc{g'^) in the case of the naive fermion with + 
Uf^ for gl = 1.2. The fine-tuned point (-2.205,0.720) as a crosspoint is located near the self- 
crossing point. 
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FIG. 10: An expanded version of Fig. [9l A blue dotted curve is the naively derived phase boundary. 
The true phase boundaries are composed of the three parts. The fine-tuned point is located slightly 
to the right and below the self-crossing point. 



We apply the same approach to the staggered Gross-Neveu model with the Adams-type 
flavored mass in Eq. ( 1251) . As seen in Eq. ( 1351) . the determinant in the logarithm in the 
effective action is given by the product of two determinants of the naive fermions with the 
mass ±Mj = ±cos(A;i/2) cos(A;2/2). Thus we only have to add the contributions from the 
two sectors. Here we take the constant a as a = — 1/a and redefine a^ + l/a — t- (Tq for a 
while. With this choice we can discuss the left cusp related to the taste with the positive 
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flavor-chirality. Then the effective potential with the (Tq shift is given by 



~ (ao-(M + l/a))^ vTo + 

•^cfflcrojTTo) = ^ 2g^ ~ ~ ' ^ 



Z,„^^^S^ + .g,.g,(-^±-^"-^)-, (69) 

, /-I ±COS%^COS^N , , 

Dt = 2ao[ ^ (70) 

We expand / with respect to Di/Dq as 

I^ = l0+J2^n, (71) 

n=l 

^^=/ T^logDj, (72) 



(n > 1), (73) 



n (27r)2 (Z}±)" 

For the continuum hmit a — 0, only the I^, if and J;^ remains nonzero as in the previous 
case. 

U + lo = Co{<yl + <) - -{<yl + vrg) log ^ (C^o = 1.177), (74) 

TT e 

/+ + /r = ^Ci (Ci = -0.896), (75) 
J2+ + I2 = -2^2 C2 (C2 = 0.404). (76) 

Details of calculations are shown in Appendix lA 2[ The effective potential and the fine- 
tuned point without 0{a) corrections {M{gl), gl{gl)) are given by the equations similar to 
Eqs. fl55|) - fl6T]) as following. The effective potential for a — in this case is given by 

'M + l/a 2 

+ f A - C'o + + - Iog4a2)ao^ + -(^o^ + tt',) log (77) 



^cff = - ( , ' + -Ci ) ao + ( — - Co + 3 log 4a2 ) TT^ 



Then the tuned point for the chiral limit without 0{a) corrections is 



20^ 

M = — ^Ci - 1, (78) 



fa 

4C2j"+l' 



9l = 777^, (79) 
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We again introduce the scale parameter (A-parameter) as 



2a A = exp 



(80) 



where we note the lattice spacing a always appears with a factor 2, which is specific to the 
staggered fermions. The coupling renormalization for the chiral and continuum limit is given 
by 

2^-^°-2^^ + ^°^(4Ab)' 
^ Co + h^g(-r^\- (82) 



2gl ^ TT ^ V4A2a2 

where we keep A finite when taking the continuum limit a — 0. Finally the renormalized 
effective potential in the chiral and continuum limit is given by 

^cfr = ^K + -o^)log^^ (83) 

In this case we take = 0.4 as an example, then the fine-tuned point is given by 

M{gl = 0.4) = -0.286, (84) 
gligl = 0.4) = 0.243. (85) 



The gap equations in this case are given by 

2 r 



(27r)2 (((Jo + ciC2)2 + vrg + s^){{ao - ciCa)^ + + 



(86) 



gl J {2nf ((ao + c.c^y + tt^ + s^){{ao - c.c^f + tt^ + s^)' 

Here we come back to the unshifted definition of o"o. In Fig. [TT]and Fig. [12] we depict the 
phase boundary M{gl) naively derived from the above gap equations for = 0.4. The 
latter is an expanded one near the self-crossing point with the true phase boundaries also 
depicted. The fine-tuned point fl8^ f l85|) is located slightly to the right and below the self- 
crossing point near the true second order phase boundary in the parity symmetric phase. 
Toward the week-coupling limit (7^ — )■ the phase structure moves down to g"^ = 0, where 
the fine-tuned point gets close to {M,g^) — (—1,0) from the parity symmetric phase. It 
means our fine-tuned point leads to the continuum theory with the chiral symmetry and 
one massless fermion corresponding to the taste with positive fiavor-chirality. The situation 
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FIG. 11: The naively derived phase boundary Mc(g^) for the staggered fermion with the Adams- 
type mass with = 0.4. The fine-tuned point (—0.286, 0.243) as a crosspoint is located near the 
self-crossing point. 

about the first order phase boundary is the same as the naive case. In Fig. [12] we depict 
the true phase boundaries for this case. In Fig. [13] we depict the order parameter ttq as a 
function of M. Here the order of the transition changes from the 2nd to the 1st around the 
order-changing point. In Fig. [Tl]we depict the (Tq potential for several values of M crossing 
the (To phase boundary. The value of ao at the minimum changes from (Tq > —1 to ctq < —1 
in a form of the Ist-order phase transition. 

We have shown that the chirally-symmetric continuum limit can be taken by fine-tuning 
a mass parameter and two coupling constants both for the naive and staggered cases. It 
indicates we obtain the two-flavor or one-flavor massless fermions in the chiral limit by 
;uning only a mass parameter when we introduce the Adams- type [l^ or Hoelbling-type 
lol flavored masses to the d = A QCD with staggered fermions. The less numerical expense 
in the staggered fermion could make the QCD simulations with these fermions faster than 
Wilson fermion. We need further investigation to answer this question. 

Now we comment on the case that we take a = in (l44l) ( |45]) . which corresponds to neither 
of the cusps but reflects the symmetries of the effective potential. At this point the coupling 
is not going to zero, and thus it is unclear how it is related to the continuum Gross-Neveu 
model. However it does seem to be possible to restore chiral symmetry there and have a 
divergent correlation length. As such it seems related to a quite special continuum limit. 
Since the M = effective potentials for the naive with M^p and the staggered fermions with 
the Adams-type mass possess the Z2 discrete chiral symmetry (ao — )■ —ctq) up to a irrelevant 
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FIG. 12: The expanded version of Fig. [TTJ A blue dotted curve is the naively derived phase 
boundary. The true phase boundaries are composed of the three parts. The fine-tuned point is 
located slightly to the right and below the self-crossing point. 
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FIG. 13: The order parameter ttq as a function of M for g'^ = 0.25, 0.26, 0.27, 0.28 where the order 
of transition changes from 1st to 2nd in Fig. [T2l 
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FIG. 14: The (Tq potential for several values of M crossing the do boundary in Fig. [T2j The value 
of do at the minimum changes from do > — 1 to (Jq < — 1 in a form of the Ist-order transition. 
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sign, the renormalization in the hnear (Tq term is prohibited. Actually we have checked Ci 
in the effective potentials as f l55|l is zero for both cases with a = 0. This is because the 
continuum chiral symmetry is broken while the discrete one is unbroken by these flavored 
masses. Thus, if we start with M = 0, it appears we need not fine-tune the mass parameter 
for the massless continuum limit with the chiral symmetry. It indicates a strange possibility 
that the chirally symmetric continuum limit of the d = A QCD with these fermions is taken 
without fine-tuning due to this symmetry. This strange situation can occur for any flavored 
mass with the discrete chiral symmetry up to a trivial sign such as Mf = cos /c^. However 
the question is whether the continuum limit stands for physically relevant theories. Indeed 
it is unlikely since the line M = is located at the same distance from the two cusps thus 
the continuum limit along it would have no physical fermions, although there might exist 
some relevant theory without fermions like the Ising theory. On the other hand, in the naive 

(2) 

fermion with Mj- or the d = A staggered fermion with the mass proposed by Hoelbling 
in the M = line has a cusp in the weak coupling region. The effective actions in 
these cases also have the discrete chiral symmetry and the same situation occurs. Thus, 
the continuum limit without fine-tuning in them may lead to the relevant theories with 
the parity symmetry being broken since the continuum limit is taken from the Aoki phase 
in these cases. This kind of the parity or CP broken theory with massless fermions would 



belong to the same universality class as minimally doubled fermions 
QCD with the sign of mass being different between the two flavors 
this topic is devoted to the future work. 



24| or the two-flavor 
Further study on 



V. SUMMARY AND DISCUSSION 



In this paper we investigate the parity-broken phase structure for naive and staggered 
fermions with the flavored mass by using the two-dimensional lattice Gross-Neveu models. 
We have shown the Aoki phase exists both in staggered and naive cases reflecting the mass 
splitting in species. 

In Sec. [Til we study the phase structure for the naive Gross-Neveu model with the flavored 
masses. We consider the two types of flavored mass terms for 2d naive fermions, which cause 
two different kinds of mass splitting in species. We also consider a linear combination of 
these terms. We solve the gap equations for the large N limit and obtain the second order 
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phase boundaries in the M-g'^ plane. The parity broken phase diagram has some common 
properties with the Wilson case, and reflects the mass splitting. We can make varieties of 
phase structures depending on arbitrary linear combinations of the two types of the masses. 
In Sec. mil we consider the generalized staggered Gross-Neveu model including two types 
of four-point interactions. We take the same process as in the case of the naive fermion 
to obtain the phase diagram for the staggered fermion with the Adams-type flavored mass. 
We show the Aoki phase exists also in this case reflecting the mass splitting of tastes. This 
elucidation can contribute to the practical application of these fermions and their overlap 
versions. In Sec. II VI we discuss the continuum limit of these Gross-Neveu models around the 
cusps in the phase digram. We show that the chirally-symmetric continuum limit with the 
number of massless species associated with each of the cusps can be taken by flne-tuning 
a mass parameter and two coupling constants in both cases. From this we speculate the 
chiral limit can be taken by flne-tuning only a mass parameter in ci = 4 lattice QCD with 
staggered fermions with the Adams-type 18] or Hoelbling-type 13] masses. It indicates 
we can obtain the one- or two-flavor massless fermions in the continuum from the staggered 
setup and regard massless pions as Goldstone bosons due to the spontaneous chiral symmetry 
breaking as in the case with Wilson fermion. These approaches avoid the use of the rooting 
approximation to reduce the number of tastes. We also study the flrst order phase boundaries 
peculiar to the two-coupling cases of the lattice Gross-Neveu models. We show there exist 
two kinds of the flrst order phase boundaries with respect to parity and chiral symmetry 
breaking as in the case of Wilson fermion. 

We comment on the possible advances of the one-flavor or two-flavor staggered fermions 
without rooting discussed in this paper compared to Wilson fermion. Taking account of 
less numerical expense in the staggered fermion, it will be numerically better than Wilson 
fermion in the lattice QCD simulations. We can estimate how good it is easily by calculating 
simple examples. Future works will be devoted to this study. 
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Appendix A: Derivation of the effective potentials 

In this appendix we evaluate the integrals which are required for the effective potentials 
for the cases with the naive and staggered fermions. 



Naive fermion 



We have to evaluate the integrals of (l50l) and (15T1) to obtain the effective potential of the 
model with the naive fermion. Let us first study the following integral. 



.la (27r)2 



log 



^0 



-1 + M/ 



(Al) 



where we denote = ^^sin^(A;^a) and Mf = cos{kia) cos(fc2a). If we omit a constant term 
which is not involving ctq and ttq, it can be rewritten in an integral representation as 



Foip) 



dp Fo(p), 



d^k 



1 



(27r)2sVa2 + (-l + M;)7a2 + p 
We pick up the divergent part in the limit of a — )■ 0, 

^3./(2a) / 1 1 

"3-/2 d^^ f 1 1 1 

-7r/2 



(A2) 
(A3) 



Fo{p) 



+ Co, 



Co 



(27r)2 \ s2 



(A4) 
0.0421). (A5) 



Here we shift the Brillouin zone to treat the divergent part, which originates from two 
massless modes around k = (0,0) and (vr,7r). We then find the following expression by 
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comparing the first term with the corresponding integral in the continuum theory. 

.3V(2a) ^ / 1 1 ^ 



— log — + Co 



0.325). (A6) 



Fo(p) = — log— +Co 

1^ ^ 



Cn = 0.367 



(AT) 



Therefore the integral is given by 

where Cq = Cq + Cg is the constant used in flS^ . By substituting this into (lA2p . we obtain 
the expression in (152]) 

1 

2^' 



/o(a ^ 0) = C{al + TTo^) - ^{al + vTo^) log ^M±Ioi 



Next we show the integral expressions of (1531) and fl54l) . They are given by 

'•^ (Pi -l + Mf 



, (27r)2s2 + (-l + M/)2 ^ 



-0.446), 



M 



f 



0.201). 



(A9) 
(AlO) 



These integrals are sufficient to consider the continuum limit of the model, but not to discuss 
the Ist-order phase transition. The 0{a) corrections come from the following integrals. 



h{a ^ 0) = = \[ 



6Ii 



a 



2(To 



-l + Mf)/ a 



-0.0923), (All) 



Mf)/a 



(27r)2 \s^/a^ + (-1 + Mff/a^ + al + tt^ s^/a^ + (-1 + Mff/a^ 



— 2(Toa 



rfp Fi(p), 



Flip) 



-1 + M;)/a 



(A12) 
(A13) 



a (2vr)2 (,Va2 + (-1 + M^) Va^ + p)^ ' 

We can evaluate the second one in a similar way by splitting into a divergent part and a 
finite constant, 



Flip) 



a-s>0 



1 r^VCSa) ^2^ 



7r/(2a) 



a) iM' 



+ 



Y.,iK-^) 



\ 



Cl 



3^2 ^2^ 



l + Mf) 



+ 



+ 



Cl 



(27r)2 W.2 + (-l + iVOT 2E,a 2E,(e.-^ 



(A14) 

0.00912). 
(A15) 
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The divergent part is given by 



^\og^ + c[ (4 = 0.166). 

(A16) 



Thus we obtain 



Fi{p) 



log- 



Ci Ci = ci 



47r "'^ a?p ' " \ " 2 
By substituting this expression into f lA12p . we obtain 



-0.0741 



5Ii = —2(TQa 



'-il^^o + ^o) + ^(^0 + ^o) log 



These integrals contribute to the 0{a) corrections to the effective potential (166|) . 



2. Staggered fermion 



(A17) 



(A18) 



We evaluate the integrals required for the effective potentials with the staggered fermion. 
Explicit forms of the finite constants in ( !75|) and ( !76|) are simply given by 



C2 




■1 + M 



-1 + M 



+ 



-l-Mf 



f) 



[-l-Mf)\ 



-1 + M 



/ 



-1 - M 



f 



-0.896), (A19) 



0.404), (A20) 



+ {-1 + Mf) J \s^ + {-l-Mf) 

where we use similar symbols as the naive fermion case, = sin2(/c^a/2), Mf = 
cos(/cia/2) cos{k2a/2). 

The integral f fMI) is slightly complicated, but can be evaluated in a similar manner. 
Omitting a constant term independent on o"o and ttq, it can be written as 



log 



- + (To + 7ro 



-1 + M 



dp F(p), 

-./a (27r)2sVa2 + (-l + M^)Va2 + p- 



(A21) 
(A22) 
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We can split this integral into a divergent part and a finite constant in the limit of a — )■ 0, 



a^o n^"" <Pk 1 



The divergent part is given by 



Thus we obtain 



" , = - log + C+' ( Co+' = 0.798 ) . (A24) 



F(p) = i log + (7o+ (^0+ = + Co+' = 0.842) . (A25) 

The corresponding integral becomes 

It {a -> 0) = Co+(ao^ + ^l) - -{al + nl) log + (A26) 

TT e 

The other integral is written as 

^ C^(al + ,rl) + 0[a), (A27) 

where we again omit a constant independent on ctq and ttq. As a result we obtain the 
expression of fl7i|) . 

/o+ + /o" = C,{al + vTo^) - 1(^0^ + 4) log + (Co = + Co" = LIT?) . (A29) 

TT e V / 
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